



Dilute Nuclear Matter in Chiral Perturbation Theory 
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We use chiral perturbation theory to compute the effective nucleon propagator in an expansion 
about low density in the chiral limit. We neglect four- nucleon interactions and focus on pion 
exchange. Evaluating the nucleon self-energy on its mass shell to leading order, we show that the 
effective nucleon mass increases by a small amount. We discuss the relevance of our results to the 
structure of compact stars. 
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As nuclear matter is compressed, eventually a transi- 
tion to quark-gluon matter occurs. By asymptotic free- 
dom, at very high densities the equation of state can be 
computed in perturbation theory [1, 2, 3]. This can be 
extended to moderate densities by various approxima- 
tion schemes [2, 3]. At low densities, the conventional 
approach is to use phenomenological potentials to fit ob- 
served properties of nuclear matter [4] , and then extrap- 
olate up in density. 

How the nuclear equation of state matches onto that 
for quark matter is of great significance for astrophysics 
[3, 5]. The standard expectation, as in Quantum Hadro- 
dynamics for example [6] , is that hadronic pressure rises 
quickly to a value near that for an ideal Fermi gas of 
quarks and even exceeds it at densities above normal 
nuclear density. In this case there is only one type of 
hadronic star, which might have a (small) quark core. 
If the hadronic pressure is small relative to that of ideal 
quarks, though, then there can be two classes of hadronic 
stars. There are "ordinary" neutron stars, which are 
mainly composed of nucleons. In addition, there are stars 
with a large quark core; their mass and radius are (ap- 
proximately) half that of ordinary neutron stars. In ad- 
dition, the pion tadpole contribution does not depend on 
the density to leading order. 

Thus it is imperative to understand the equation of 
state for nuclear matter. In Ref. [7], Savage and Wise 
compute mass shifts for the baryon octet using chiral 
perturbation theory, including all operators which con- 
tribute to leading order in the density. Due to four- 
nucleon interactions, they find that all masses decrease 
with increasing density; extrapolating to nuclear matter 
densities, the shifts are considerable. The self-energies 
were computed at zero momentum, though, while the 
physical point is on the mass shell. In the presence of 
a Fermi sea, the mass shell changes, in a way which is 
easily computed. To leading order in the density, the 
difference in mass shell only affects exchange, and not 



contact, terms. In this paper we compute the nucleon 
self energy on its mass shell, from the diagram for pion 
exchange. In this case, unlike Ref. [7], the usual loga- 
rithms of chiral perturbation theory appear on the mass 
shell. We note that the shift in the nucleon mass is dom- 
inated by contact terms, not single pion exchange, so the 
following exercise is a minor point of principle. 

We assume that nucleons, of mass m, are heavy, and 
that pions, with mass m n , are very light. The interaction 
of pions and nucleons is determined by the spontaneous 
breaking of chiral symmetry [6, 8, 9, 10, 11]. For light 
pions, to leading order in chiral perturbation theory, the 
only parameter which enters is the pion decay constant, 
U w 93 MeV. 

The proceeding calculation is elementary, and, besides 
those of Savage and Wise [7], it is similar to computa- 
tions by Horowitz and Serot [9], by Bernard, Kaiser, and 
Meifiner [10], and by Meifiner, Oiler, and Wirzba [11]. 

The parameters which enter can be understood with- 
out explicit computation. A Fermi gas of nucleons is 
characterized by a Fermi momentum, pf\ up to and in- 
cluding nuclear matter densities, pf <C m. In an expan- 
sion about low densities, the natural parameter which 
enters into the nucleon propagator is just the density, 
Tinucl ~ Pf- We would like a dimensionless parameter 
to characterize the expansion. At one loop order, chi- 
ral perturbation theory brings in two powers of l/f w - 
The only other parameter in the problem is the nucleon 
mass, m (at least for vanishing pion mass). Thus to lead- 
ing order, the corrections to the nucleon propagator are 
proportional to 



(1) 



which we now compute. 

To leading order in chiral perturbation theory, we take 
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the nudeon Lagrangian to be 



7 n + m ■ 



■9A n. 



(2) 



where fj, = ^Jp 2 + m 2 is the chemical potential, qa ~ 1.2 

is the axial vector coupling constant, and 7r = 7r a cr a /2, 
where the cr°'s are Pauli matrices in SU(2) flavor. Other 
interactions, such as between two nucleons and more than 
two pions, involve more powers of and so enter be- 

yond leading order in the density. There are tadpole con- 
tributions, with a pion in the loop, but like the tadpoles 
from four-nucleon interactions, these arc independent of 
the external momentum, and so of the choice of mass 
shell. 

Thus we consider single pion exchange, which con- 
tributes to the nucleon self energy, S, as 



E(P) = -: 



3g\ f d 4 K 



Af 2 J (2tt)4 {P-Kf+m 2 



75 (/>- fC) 
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—i]K — 7°^i + m 



75 (/>- JO; 
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P = (p,p) is the four-momentum of the nucleon. The 
diagrams are evaluated using the imaginary time formal- 
ism. 

In general, the nucleon self energy is a rather 

complicated function of p° and p [9, 10]. Here we shall 
only compute the nucleon self energy at a special point, 
on its mass shell: 



Pms = i (p- E p) 



(p 2 



f 



P 2 ) 
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E p is the energy of a nucleon with momentum p, so [i = 
E Pf , and E p « m + p 2 / (2m) + . . .. We also assume that 
the momentum p is on the order of the Fermi momentum, 
but it need not be especially near p f . 

Working on the mass shell allows us to greatly simplify 
the calculation. As we are working at nonzero fermion 
density, it is convenient to do the integral over k° first, 
and then integrate over k. In the imaginary time formal- 
ism, we first compute the diagram for real p°, and then 
analytically continue to imaginary values of p°, as in (4). 

In the integrand, there are four poles: two from the 

pion propagator, at k° = p° ±i\J (k - p) 2 + m 2 , and two 
from the nucleon propagator, at k = ± Ek)- Closing 
the contour in imaginary p° plane, only those poles in 
the upper half plane contribute. 

All we are interested in, though, are the density depen- 
dent effects. Of the four poles in the one loop diagram for 
the nucleon propagator, clearly one is special. The pole 
at which fco = i {p — Ek) is in the upper half plane when 
k < pf, and moves into the lower half plane when k > Pf- 
For the other three poles, the sign of their imaginary part 
does not change with k. 



All of the density dependent effects in the one loop 
diagram for the nucleon propagator are due to the shift 
in this one pole. To see this, note that the chemical 



potential [i only enters by changing ip° 



ip" 



/i. If we 



work on the mass shell, however, then ip° + fi = +E p : 
while p° changes with fi, ip° + [i does not. Thus if no 
poles switched the sign of their imaginary part, then we 
would find that there were no density dependent effects 
in the propagator at one loop order. For example, there is 
wave-function renormalization for the nucleon field, but 
given that [i enters just as a shift in p°, this is standard; 
there is no new wave- function renormalization associated 
with ^7°, separate from p 1 . 

The contribution of the pole at fc° = i (p — Ek) is sim- 
ple to include: one only integrates over momentum below 
the Fermi surface. A similar result is found, rather more 
immediately, using the real time formalism. There, the 
nucleon propagator is the sum of two terms, one the same 
as in the vacuum, plus a density dependent term. 

To pick up the contribution of just this one pole, we 
take 



dk° 



1 



2tt -i JC - 7>- 



m 



2E, 



(5) 

On the right hand side, k n = ifai — Ek), and only |fc| < pf 

contribute to the integral over k. 

Now we need to sandwich the inverse nucleon propa- 
gator in this expression between the 75 {p— fi)'s from 
the pion vertices. There are three types of terms which 
contribute. One is from the term ~ m in the nucleon 
propagator: 

cm - 75 ( P- ^)m 75 ( P- fC) = -m(P - K) 2 , (6) 
one from the term ~ /17 : 

c M = 75 (f- ^)M7°75(/'- K) (7) 

= 7 V ((P° - fc°) 2 - (P ~ kf) + 2 M (p° - k°) (f- f) , 
and one from the term ~ /T: 

c p = 75 (P- fC) i ^75 (P- fC) (8) 

= i (2(P -K)-K P-(P 2 - K 2 ) fC) . 

To compute the leading terms about small density, we 
can greatly simplify these expressions. For example, for 
the nucleon energy, we can replace Ek ~ m, since correc- 
tions are down by (pf/m) 2 . Further, as we are comput- 
ing on the mass shell, the energy p° is small relative to 
the spatial momentum; in magnitude, as p° ~ p 2 /m, 
p° is down by Pf/m relative to p. This means that 
in the pion propagator, and in c m , (6), we can replace 
(P-K) 2 w {p-k) 2 . 

For the other contributions, one must be careful to 
keep track of relatively small terms, and also ~ p 7 . 
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For c M , (7), for the piece 7 we can drop (p° — fc ) 2 relative 
to (p— fc) 2 , and take « m. However, for the piece 
~ n(p° — fc°), we have to keep track of the subdominant 
term, so 

c M » - m 7 ° (p - fc) 2 + i (p 2 - fc 2 ) I) . (9) 

For the last term, c p in (8), we can approximate 

c p « i (2(p - fc) • fc f - (p 2 - fc 2 ) jr) . (10) 

We keep the terms ~ 7 , which are nominally down by 
Pf/m, in order to extract the term ~ p°7°. 

Adding all of these terms together, we find a remark- 
able simplification: 

S(pL.P) « - ((iPm. + M) 7° + * f+ m) S (p), (H) 
where 

51 , , 3^ y ^fc / ( P --fc) 2 \ 

We also checked that the same result is found using the 
real time formalism. 

This form is illuminating, because it is obvious that 
in the chiral limit, when m„ = 0, the function S is 
independent of momentum: 

So(rt - + i6^wf- + ^^F' (13) 

where n nuc i = 2p^/(37r 2 ) is the density of nucleons. 

Away from the chiral limit, ^ 0, So is momentum 
dependent: 

s *>=Ti^4( l+ W ^s ° <,,, )' <14) 

V 4 PP/ y v(p/+p) 2 + "»jy 

+ ^ farctan f + arctan f EtZ£^ - 1 

(15) 

At zero momentum, p = 0, this agrees with Savage and 
Wise [7]. We see that chiral logarithms appear when 
p^O, although there is a arctan(p//m) at p = 0. These 
chiral logarithms are standard [8] , and relatively innocu- 
ous. Even at the Fermi surface, p = Pf, they vanish like 
to 2 log(m 7T ) as m^r — > 0. 

One can compute So as a function of m,. For illus- 
tration, consider its value at the Fermi surface. Then 
one can show that increasing the pion mass tends to de- 
crease the value of So; as m n — > 00, Sq vanishes like 



« p^/(m/ 2 )(p 2 /m 2 ). That S vanishes like w 1/m 2 at 
large m T is evident from the integral representation, (12). 

We can use these results to compute the nature of nu- 
cleon quasiparticles. Adding the self energy, the effective 
nucleon propagator is 

KffiPLiP) = ^bare £ (16) 

= - {{iP° ms + M) 7° + if) (1 - So) + m (1 + S ) . 

The change in the position of the pole in the nucleon 
propagator is easy to compute. In particular, the mass 
of the nucleon is shifted up: 

me//=m (rr§) »"»(1 + 2S ) • (17) 

This expression holds in the chiral limit. Comparing with 
(9), half of the mass shift arises from the shift in the term 
~ m, and half from what can be viewed as wave-function 
renormalization. 

Away from the chiral limit, where So is a function of 
momentum, the change in the mass cannot be read off so 
immediately. In that case, one has to define the effective 
mass by other means, as in (11.66) of [13]. 

This increase in the effective nucleon mass is in con- 
trast to what happens at zero density, but nonzero tem- 
perature. To leading order in an expansion about zero 
temperature, in the chiral limit the nucleon mass does 
not shift to - T 2 [14]. 

At normal nuclear matter density, pf « 270 MeV. 
The correction which we computed, from single pion ex- 
change, is tiny, 2S ~ .04. This suggests that chiral 
perturbation theory might be a reasonable guide to the 
properties of nucleons, even at nuclear matter densities. 

This conclusion is premature. While the corrections to 
the nucleon propagator are very small, corrections to the 
pion propagator are large. For most momentum, such 
as near the pion mass shell, the corrections to the pion 
propagator are like those of the nucleon, proportional to 
the density, ~ p^/(m/ 2 ), (1). If the pion is far off its 
mass shell, though, with an energy ui ~ p 2 /m, it can 
scatter into a nucleon particle-hole pair. For such nearly 
static pions, the pion self energy is enhanced by a factor 
of m/u) <~ m 2 /p 2 . The correct expansion parameter for 
the pion propagator is then not Pj/(to/ 2 ), but 

Pf ™ 2 ( 9a \ m Pf M8 x 

^ri ~% ~ U^J ir • (18) 

Numerically, this parameter is much larger than S in 
(13). In fact, as we are dealing with a non-relativistic 
system, this enhancement of the pion propagator is well 
known from condensed matter physics, and represents 
the need to rcsum the nearly static pion propagator 
through the Random Phase Approximation (RPA) [13]. 
Indeed, the factor of g^/(27r 2 ) arises from an explicit cal- 
culation in the RPA limit, from (4.21) of Meifiner, Oiler, 
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and Wirzba [11]. For normal nuclear matter density, the 
parameter of the RPA pion propagator in (18) is w 2 
at nuclear matter densities. Since this parameter is only 
linear in the Fermi momentum, if we require that this 
parameter be less than, say, 1/2, this means that we can 
use chiral perturbation theory to compute the nuclear 
equation of state only up to p/ ~ 70 MeV. This cor- 
responds to densities which are (1/4) 3 = 1/64 those of 
normal nuclear matter! 

This restriction on the use of chiral perturbation theory 
is not that surprising. In computing the free energy, the 
typical pion momentum is of order ~ pf, with energies 
~ Pf/m. To use a chiral Lagrangian, the pion momen- 
tum should be small relative to f v , which is similar to 
the condition derived from (18). What is not evident is 
while there is a factor of l/(27r 2 ) from chiral perturbation 
theory in (18), this is compensated by the factor of the 
nucleon mass in the numerator. 

Nevertheless, such computations [10, 11] are manifestly 
of interest, so as to gain a more general understanding of 
the nuclear equation of state. Carrying out such calcula- 
tions beyond leading order is technically very challenging. 
Using an RPA corrected propagator for the nearly static 



pion is straightforward. What is difficult is knowing how 
to separate diagrams with two pion exchange from other 
effects. In Quantum Hadrodynamics [6], one must sepa- 
rate two pion exchange from that of of heavier mesons, 
such as the a and the u>. In "pionlcss" effective theories 
[15], two pion exchange contributes to point like interac- 
tions between four or more nucleons. 

We can draw some tentative conclusions about the 
hadronic pressure, which motivated this study. To lead- 
ing order, the non-ideal terms in the pressure are propor- 
tional to S , which is very small [10]. At higher order, 
even if corrections to the pion propagator are large, their 
effect on the nucleon propagator, and the free energy, 
can still be small, as a large correction to a small num- 
ber. Thus the possibility of a hadronic phase with a small 
pressure, required for a new class of quark stars, remains 
viable. 

Acknowledgements: We thank E. Kolomeitsev, M. 
Savage, and U. van Kolck for discussions. E.S.F. is 
partially supported by CAPES, CNPq, FAPERJ and 
FUJB/UFRJ. J.S.B. is partially supported by DOE grant 
DE-FG-02-93ER-40764. The research of R.D.P. is sup- 
ported by DOE grant DE-AC02-98CH10886. 



[1] B. A. Frccdman and L. D. McLerran, Phys. Rev. D16, 
1130 (1977); ibid., D16, 1147 (1977); ibid., D16, 1169 
(1977); ibid., D17, 1109 (1978); V. Baluni, Phys. Rev. D 
17, 2092 (1978). 

[2] R. Baier and K. Redlich, Phys. Rev. Lett. 84, 2100 
(2000); hep-ph/9908372; J. P. Blaizot, E. Iancu and 

A. Rebhan, Phys. Rev. D 63, 065003 (2001), hep- 
ph/0005003; A. Peshier, B. Kampfer and G. Soff, Phys. 
Rev. D 66, 094003 (2002), hep-ph/0206229. 

[3] E. S. Fraga, R. D. Pisarski, and J. Schaffner-Bielich, 
Phys. Rev. D 63, 121702 (2001), hep-ph/0101143; Nucl. 
Phys. A 702, 217 (2002), nucl-th/0110077. 

[4] A. Akmal, V. R. Pandharipande and D. G. Ravenhall, 
Phys. Rev. C 58, 1804 (1998), nucl-th/9804027, and ref- 
erences therein. 

[5] U. H. Gerlach, Phys. Rev. 172, 1325 (1968); N. K. Glen- 
denning and C. Kettner, Astron. Astrophys. 353, L9 
(2000), astro-ph/9807155; D. Blaschke, H. Grigorian, 
G. Poghosyan, C. D. Roberts and S. M. Schmidt, Phys. 
Lett. B450, 207 (1999), nucl-th/9801060; A. Peshier, 

B. Kampfer and G. Soff, Phys. Rev. C61, 045203 (2000), 
hep-ph/9911474. 

[6] B. D. Serot and J. D. Walecka, Intl. Jour. Mod. Phys. 
6, 515 (1997), micl-th/9701058; R. J. Furnstahl and 
B. D. Serot, Comments Nucl. Part. Phys. 2, A23 (2000), 
nucl-th/0005072. 

[7] M. J. Savage and M. B. Wise, Phys. Rev. D 53, 349 
(1996), hep-ph/9507288. 

[8] J. F. Donoghue, E. Golowich, and B. R. Holstein, "Dy- 
namics of the Standard Model" (Cambridge University 



Press, Cambridge, 1992). 
[9] C. J. Horowitz and B. D. Serot, Phys. Lett. B 108, 377 
(1982). 

[10] V. Bernard, N. Kaiser, Ulf-G. Meifiner, Intl. Jour. 
Mod. Phys. E4, 193 (1995), hep-ph/9501384; M. Lutz, 
B. Friman, and C. Appel, Phys. Lett. B 474, 7 (2000), 
nucl-th/9907078; N. Kaiser, S. Fritsch, and W. Weise, 
Nucl. Phys. A 697, 255 (2002), nucl-th/0105057; Nucl. 
Phys. A 700, 343 (2002), nucl-th/0108010; Phys. Lett. B 
545, 73 (2002), nucl-th/0202005; nucl-th/0212049. 

[11] Ulf-G. Meifiner, J. A. Oiler, and A. Wirzba, Annals Phys. 
297, 27 (2002), micl-th/0109026. 

[12] G. E. Brown and M. Rho, Phys. Rev. Lett. 66, 2720 
(1991); Phys. Rep. 363, 85 (2002), hep-ph/0103102. 

[13] A. L. Fetter and J. D. Walecka, "Quantum Theory 
of Many-Particle Systems" (New York, McGraw-Hill, 
1971). 

[14] H. Leutwyler and A. V. Smilga, Nucl. Phys. B 342, 302 
(1990). 

[15] S. Weinberg, Nucl. Phys. B 363, 3 (1991); D. B. Ka- 
plan, M. J. Savage and M. B. Wise, Nucl. Phys. B 478, 
629 (1996), nucl-th/9605002; ibid. 534, 329 (1998); nucl- 
th/9802075, U. van Kolck, Prog. Part. Nucl. Phys. 43, 
337 (1999), nucl-th/9902015, J. W. Chen, G. Rupak and 
M. J. Savage, Nucl. Phys. A 653, 386 (1999), nucl- 
th/9902056, S. R. Beane, P. F. Bedaque, M. J. Savage 
and U. van Kolck, Nucl. Phys. A 700, 377 (2002), nucl- 
th/0104030. 



